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Abstract 

Mass/electron  transports  in  a  passive  cathode  of  a  proton  exchange  membrane  (PEM)  fuel  cell  have  been  studied  numerically.  The  porous  cathode 
in  contact  with  a  perforated  current  collector  breathes  fresh  air  through  an  array  of  orifices.  Diffusions  of  reactant  species  in  the  porous  cathodes 
are  described  by  the  Stefan-Maxwell  equation.  Electrochemical  reaction  on  the  surfaces  of  the  porous  cathode  is  depicted  via  the  Butler- Volmer 
equation.  Gas  flow  in  the  air-breathing  porous  cathodes  is  governed  by  isotropic  linear  resistance  model  with  constant  porosity  and  permeability. 
The  electron/ion  transports  in  the  catalyst/electrolyte  are  dealt  with  the  charge  conservations  based  on  the  Ohm’s  law.  A  finite-element  method  is 
employed  to  solve  the  above-coupled  equations.  The  effect  of  overpotential  on  the  fluid  flow,  mass  transport  and  electrochemistry  is  examined. 
Detailed  electrochemical/mass  characteristics  such  as  flow  velocities,  species  mass  fraction,  species  flux  and  current  density  distributions  are 
presented.  They  can  provide  a  solid  basis  for  optimizing  the  geometry  of  the  PEM  fuel  cell  stack  running  with  a  passive  mode. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  increasing  number  of  portable  electronic  devices  on  the 
market  today  such  as  laptop  computers,  cellular  phones,  PDAs, 
remote-controlled  toys,  is  creating  a  demand  for  improved,  more 
environmentally  friendly  technologies  in  energy  storage  and 
conversion.  Small-scale  fuel  cells  are  regarded  as  a  potential 
candidate  in  replacing  batteries  in  portable  applications  [1,2]. 
However,  up  to  the  present  time,  fuel  cells  were  still  too  com¬ 
plicated  and  expensive  to  meet  the  need  because  they  required 
cooling,  humidification  and  pressurization  sub-systems  in  opera¬ 
tion.  In  response,  researchers  and  engineers  are  trying  to  develop 
passive  fuel  cells  for  portable  power  applications.  One  interest¬ 
ing  approach  is  air-breathing  fuel  cells,  i.e.  the  oxygen  needed 
by  the  electrochemical  reaction  is  taken  directly  from  the  sur¬ 
rounding  air  by  diffusion  and/or  natural  convection.  The  fuel  cell 
running  in  a  passive  manner  has  the  benefits  of  minimizing  the 
number  of  power  consumed  auxiliary  devices  that  enhances  the 
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system  efficiency.  However,  the  power  density  of  cell  itself  is  not 
as  good  as  that  running  in  an  active  mode  due  to  poor  transport 
mechanisms  [3-13].  It  is  well  known  that  the  current  production 
by  electrochemical  reaction  is  directly  proportional  to  the  local 
oxygen  concentration  in  the  fuel  cell.  Inadequate  airflow  cannot 
provide  enough  oxygen  for  the  electrochemical  reaction  on  the 
active  surfaces.  It  results  in  heterogeneous  current  distribution 
in  the  electrode  that  reduces  the  performance  of  the  fuel  cell. 
Thus,  one  of  the  most  challenges  in  the  design  of  passive  fuel 
cells  is  how  to  feed  the  electrode  with  reactants  and  remove  the 
products  from  the  electrode  appropriately. 

The  objective  of  the  present  paper  is  to  numerically  study 
the  transport  phenomena  in  a  passive  mode  cathode  of  a  PEM 
fuel  cell  that  operates  without  external  forced  flows.  The  porous 
cathode  is  attached  to  a  perforated  current  collector.  It  breathes 
the  fresh  air  through  the  perforations  and  routes  out  the  electric 
current  through  its  solid  counterpart.  A  set  of  conservation  equa¬ 
tions  of  mass,  momentum,  species  and  charge  is  developed.  The 
Darcy  law  describes  the  fluid  flow  characteristics  in  the  porous 
cathode.  The  Stefan-Maxwell  equations  illustrate  the  multi¬ 
species  diffusion  in  the  porous  cathode.  The  Butler- Volmer 
equation  depicts  the  electrochemical  reaction  on  the  surfaces  of 
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Nomenclature 


Ci 


Dh2o 

Dq2 

F 

Ft 

M 

n 

P 

Pi,m 


Pi,s 


R 

T 

u 


x,  y,  z 


mole  concentration  of  the  species  i  (molm  3), 

N 

r-  —  2/L  f\  '  0Ji  P 

Ll  ~  Mi'  Z-^Mj  RT 

j 

diffusivity  of  water  vapor  (m2  s-1) 
diffusivity  of  oxygen  (m2  s-1) 

Faraday’s  constant 
current  density  (A  m-2) 
molecular  weight  (kg  mol-1) 
mass  flux  (kgm2  s-1) 
pressure  (Pa) 

possibility  of  the  electrolyte  in  the  connection  of 
the  catalyst  layer 

possibility  of  the  catalyst  in  the  connection  of  the 
catalyst  layer 

universal  gas  constant  (W  mol-1  K-1) 
temperature  (K) 
velocity  vectors  (ms-1) 
coordinate  system,  Fig.  1  (m) 


Greek  symbols 
a  symmetric  factor 

s  cathode  porosity 

<p\  potential  of  the  ionic  conductor  (electrolyte- 

phase)  (V) 

0s  potential  of  the  electric  conductor  (catalyst- 

phase)  (V) 

rjt  total  overpotential  across  the  computational  mod¬ 
ule  (V) 

k  permeability  (m2) 

li  dynamic  viscosity  (m  s-2) 

p  density  (kgm3) 

g\  ionic  conductivity  of  the  membrane  phase 

(S  m-1) 

crs  electric  conductivity  of  the  catalyst-phase 

(S  m-1) 

r  tortuosity 

v\  volume  fraction  of  the  ionic  conductor 

(electrolyte-phase)  in  the  cathode 
vs  volume  fraction  of  the  electronic  conductor 

(catalyst-phase)  in  the  cathode 
co  mass  fraction 


the  porous  cathode.  A  finite-element  based  computational  fluid 
dynamics  (CFD)  [14,15]  methodology  is  employed  to  solve  the 
multi-physics  transports  in  the  air-breathing  cathode  of  a  PEM 
fuel  cell.  Emphasis  is  placed  on  the  effect  of  the  total  overpoten¬ 
tial  on  the  species/electron  co-transports  in  a  porous  cathode  of 
a  PEM  fuel  cell.  The  relation  between  the  current  distributions 
and  the  mass  transports  are  discussed  in  detail.  The  informative 
results  obtained  by  the  present  study  can  help  in  understanding 
of  the  local  gas  transports  and  electrochemical  characteristics  in 
an  air-breathing  PEM  fuel  cell.  Also,  they  can  provide  a  solid 


Reactive  layer  - 
(Porous  cathode) 
Electrolyte  layer 
(PEM)  ' 


Computational  module 


Breathing  holes 
(Gas  inlet  and  outlet) 


Fig.  1.  Fuel  cell  cathode  with  a  perforated  current  collector. 


basis  for  optimizing  the  geometry  of  the  passive  PEM  fuel  cell 
stack. 

2.  Model  descriptions 

2.7.  Governing  equations 

Fig.  1  shows  a  typical  geometry  of  a  self-breathing  cathode 
of  a  PEM  fuel  cell  with  perforated  current  collectors,  which  is 
often  seen  in  mini-  or  micro-fuel  cells  [16,17].  The  pink  circular 
zones  on  top  surfaces  of  the  module  are  the  cathode  inlets  that 
are  open  to  the  ambient,  while  the  rest  of  the  top  surfaces  sit 
flush  against  the  current  collector.  Due  to  the  perforation  layout, 
a  3D  model  is  needed  to  simulate  the  mass  transport,  current  and 
reaction  distributions.  For  memory  and  time  saving,  a  unit  cell 
shown  in  Fig.  2  is  considered  in  the  present  study.  The  unit  cell 
is  1.5  mm  x  1.5  mm  in  surface,  and  the  gas  inlet/outlet  hole  has 
a  radius  of  1.0  mm.  Both  the  reactive  layer  and  the  electrolyte 
layer  are  75  pm  thick.  The  origin  of  the  coordinate  system  is  set 
to  the  remote  comer  of  the  unit  cell,  as  shown  in  Fig.  2. 

The  upper  rectangular  domain  is  a  porous  structure  that  con¬ 
tains  the  feed  gas  mixture,  an  electronically  conducting  material 
(catalyst)  and  an  ionically  conducting  material  (electrolyte).  The 
oxygen  reduction  reaction  takes  place  in  the  reactive  layer.  It  can 
be  represented  by  the  following  equations. 

02  +  4H+  +  4e“  —*■  2H20  (1) 
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Note  that  the  circular  holes  in  the  collector  are  not  only  the 
entrance  of  the  fresh  air  but  also  the  exit  of  the  products  (water 
vapor)  of  the  electrochemical  reaction. 

The  lower  domain  corresponds  to  an  electrolyte  layer  (PEM), 
interconnecting  the  two  electrodes  of  the  fuel  cell.  No  reaction 
takes  place  in  this  layer.  Neither  are  there  pores  to  allow  gas 
to  flow  nor  material  for  electronic  current.  The  ionic  current 
balance  in  the  electrolyte  membrane  between  the  anode  and  the 
cathode  is  given  by: 

V(-oi  V0O  =  0  (2) 

where  a  \  is  the  ionic  conductivity  of  the  electrolyte  and  0 1  is  the 
electrolyte  potential. 

Two  kinds  of  solid  materials  consist  of  the  porous  matrices 
of  the  cathode  i.e.,  electrolyte  (ionic  conductor)  and  catalyst 
(electronic  conductor).  A  potential  difference  exists  between  the 
electrolyte  and  the  catalyst  to  drive  the  transfer  current,  keeping 
the  electrochemical  reaction  continuously.  The  current  passes 
through  cathode  can  be  decomposed  two  parts,  i.e., 


V  •  (pU6t»H2o)  =  pT>H20,eff  Vu;o2  + 


pco  o2 

M 


Dh2o  VM  +  Su2o 


(11) 


<^N2  =  1  —  <J>02  —  OJVh_0 


(12) 


where  p  is  the  density  of  the  gas  mixture,  u  the  velocity  vector 
and  M  is  the  mole  masses  of  the  mixture,.  The  first  term  on 
the  right-hand  side  of  Eqs.  (10)  and  (11)  represents  the  Fickian 
diffusion  due  to  concentration  gradient.  While  the  second  one 
on  the  right-hand  side  of  Eqs.  (10)  and  (1 1)  is  the  correction  to 
enforce  the  Stefan-Maxwell  equations  for  the  multi-component 
diffusion  system.  The  source  terms  of  the  mass  balances  for 
oxygen  and  water  vapor  are,  respectively,  given  by: 


At5aMo2 
4 F 


(13) 


/ct*Sa^H2Q 

IF 


(14) 


i  =  h  +  h  (3) 

4  and  i\  are  the  currents  flowing  through  the  catalyst  and  the 
electrolyte,  respectively.  Since  the  cathodes  are  electroneutral 
everywhere,  there  is  no  charge-buildup  in  the  catalyst  layers. 
Thus,  the  charge  conservation  is 

V  •  i  =  0  (4) 

That  is 

V  ■  is  =  -  V  •  h  (5) 

These  two  current  components  interact  through  electrochemical 
reactions.  The  electrons  are  transferred  to  the  electrolyte  from 
the  catalyst.  Application  of  Ohm’s  law  to  equation  yields  the 
current  conservation: 

V  ’  (— °s,eff  V0S)  —  —  ’  At  (6) 

^  ’  (— ^Teff  V0i)  =  ’  At  (2) 

where  0S  and  0 1  are  the  catalyst  potential  and  the  electrolyte 
potential,  respectively.  Sa  is  the  surface  area-to-volume  ratio, 
and  /ct  is  the  local  transfer  current  density.  <rs,eff  and  a  peff  are 
the  effective  electronic  and  ionic  conductivities  of  the  catalyst 
and  electrolyte,  respectively.  They  are  modeled  as 


where  F  is  the  Faraday’s  constant.  The  above  effective  diffu- 
sivities  follow  the  Bruggemann  model  [13]  that  describes  the 
species  diffuses  in  the  porous  electrode  of  a  multi-component 
mixture  system,  i.e., 

Ao2,eff  =  stDo2  (15) 

£>H20,eff  =  STDn20  (16) 


The  velocity  vector  is  coupled  to  Darcy’s  law  accordingly: 
u  =  --Wp  (17) 

where  k  is  the  permeability,  (i  the  viscosity  and p  is  the  pressure. 
The  continuity  equation  for  the  gas  flow  mixture  is  the  sum  of 
the  continuity  for  the  three  involved  species,  which  yields: 


V 


P-Vp) 

A i  J 


At  5a 47  At5a47H20 

4  F  IF 


(18) 


The  local  charge  transfer  current  density  couples  the  current 
balances  with  the  mass  balances,  according  to: 


co2  \ 

-  exp 

c02,ref  / 


4a  F 
~RT 


OAeff  —  ^(l  Sc)  X  Us  X  pi  s 
01,eff  =  Ol(l  -  Sc)  X  V\  X  Pi  i 


(8) 

(9) 


^  cu2o 


CH20,ref 


2 

exp 


4(1  —  a)F 
RT 


(19) 


where  vs  and  v\  are  the  volume  fraction  of  the  catalyst  and 
electrolyte  in  the  catalyst  layer,  respectively.  p^s  and  p^\  are  the 
possibilities  of  the  catalyst  and  electrolyte  in  the  connection  of 
the  catalyst  layer,  respectively  [18,19].  It  is  noted  that  only  a 
long-range  connection  of  the  same  particles  stretch  from  the 
catalyst  layer  to  the  electrolyte  ensures  good  conductivity. 

The  conservations  of  species  for  oxygen,  water  vapor  and 
nitrogen  in  the  gas  phase  are,  respectively,  expressed  as 


V  •  ( p\xcoo2 )  =  pDo2  eff  V&>o2  + 


poi  o2 


Dq2  VM  +  Sq, 


(10) 


where  a  is  the  symmetric  factor,  R  the  universal  gas  constant 
and  T  is  the  temperature.  co2  and  ch2o  are  the  concentrations  of 
oxygen  and  water  vapor,  respectively. 

2.2.  Boundary  conditions 

The  boundary  conditions  for  the  ionic  current  balances  are 
insulating  at  all  boundaries,  i.e., 

n  •  (— oAeff  V0i)  =  0 


M 


(20) 
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except  for  the  lowest  boundary  for  the  electrolyte  layer.  The 
condition  at  this  boundary  sets  the  ionic  potential: 

01  =  0i,o  (21) 

The  boundaries  for  the  electronic  current  balance  are  those  that 
encompass  the  cathode  sub-domain.  The  boundary  conditions 
are  all  insulating  except  for  the  current  collector  boundary  where 
the  potential  is  set  to  zero  (ground  connection): 

0s  =  0  (22) 

The  boundary  conditions  for  the  species  balances  are  also  those 
that  encompass  the  electrode  sub-domain  and  are  all  insulating 
except  for  the  gas  channel  boundary.  The  conditions  at  the  mod¬ 
ule  inlet  (perforations)  are  given  by  the  composition  of  the  inlet 
gas 

o>i  =  o>u  0  (23) 

For  the  momentum  balance  through  Darcy’s  law,  all  boundaries 
are  insulating: 

n-fSvpj=°  (24) 

except  for  the  module  inlet,  where  the  pressure  is  set: 

P  —  Patm  (25) 

The  input  data  to  the  model  is  given  in  Table  1 . 

2.3.  Numerical  schemes 

The  governing  equations  are  numerically  solved  by  the  finite- 
element  method  [20-22].  It  uses  the  Broy den’s  method  with 
an  LU-decomposition  pre-conditioner  to  solve  the  non-linear 
equations  iteratively.  To  reduce  continuity  errors,  a  penalty  term 
is  employed  for  pressure.  Thus,  there  is  a  continuous  part  of 
the  pressure  and  piecewise  constant  part  providing  and  extra 
degree  of  freedom  (DOF)  for  pressure  on  each  element.  It  uses 
Newton-Raphson  iteration  to  solve  the  close-coupled  groups 

Table  1 


Electrochemical  and  flow  characteristic  data  used  in  the  present  model 


Expression 

Symbols 

Data 

Exchange  current  density  (Am2) 

io 

1.0  x  icr3 

Surface  to  volume  ratio  (m-1) 

sa 

1.0  x  109 

Ambient  temperature  (K) 

T 

353 

Electric  conductivity  of  the  catalyst  (S  m- 1 ) 

1000 

Ionic  conductivity  of  the  electrolyte  (S  m- 1 ) 

CTl 

5.01 

Dynamic  viscosity  of  gas  mixture  (ms-2) 

/X 

3.0  X  IO"5 

Porosity  of  the  cathode 

e 

0.4 

Permeability  of  the  cathode  (m2) 

K 

1.0  X  IO"9 

Reference  mole  fraction  of  oxygen 

c02,ref 

3.6641 

Reference  mole  fraction  of  water 

c02,ref 

3.6641 

Ambient  pressure  (Pa) 

Pin 

1.013  x  105 

N2  molecular  weight  (kg  mol-1) 

MN2 

0.028 

O2  molecular  weight  (kg  mol- 1 ) 

mo2 

0.0321 

H2O  molecular  weight  (kg  mol-1) 

Mh2o 

0.018 

Mass  fraction  of  O2  at  inlet 

(^,02,in 

0.233 

Mass  fraction  of  N2  at  inlet 

&>N2,in 

0.766 

Mass  fraction  of  H2O  at  inlet 

^H2  O,  in 

0.001 

(velocity,  pressure,  temperature,  concentration  and  electricity) 
and  uses  the  frontal  algorithm  (Gaussian  elimination)  to  solve 
the  linearized  system  of  equations  for  each  iteration. 

In  the  computational  domain,  a  total  of  1828  nodes  and  8789 
meshes  were  used  (quadratic  velocities  in  each  direction).  Addi¬ 
tional  runs  for  the  coarser  meshes,  7712,  and  the  finer  meshes, 
10982,  are  taken  for  a  check  of  grid  independence.  The  pro¬ 
gram  gives  results  within  1  %  of  each  other  on  the  finest  meshes 
used.  A  comparison  of  the  results  of  the  two  mesh  sizes,  8789 
and  10982,  shows  that  the  maximum  discrepancies  in  the  axial 
velocity  and  oxygen  concentration  profiles  are  1.5  and  1.6%, 
respectively.  In  addition,  results  indicate  a  maximum  change 
of  2.1%  in  current  density  distribution  between  the  solutions 
of  8789  and  10982  meshes.  These  changes  are  so  small  that  the 
accuracy  of  the  solutions  on  8789  meshes  is  deemed  satisfactory. 
Solutions  are  considered  to  be  converged  at  each  test  condition 
after  the  ratio  of  residual  source  (including  mass,  momentum  and 
species)  to  the  maximum  flux  across  a  control  surface  becomes 
below  1.0  x  10“ 6.  The  CPU  time  ranged  from  100  to  500  min 


(c) 

Fig.  3.  Effect  of  overpotential  on  the  velocity  distributions:  (a)  rjt  =  0.35  V,  (b) 
rjt  =  0.55V  and  (c)  rjt  =  0.70  V  (unit:  ms-1). 
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Fig.  4.  Iso-surface  distributions  of  oxygen  mass  fractions  and  water  vapor  mass 
fractions,  r) t  =  0.55  V. 


(mm) 


(mm) 


Diffusive  mass  flux  of  oxygen 


Convective  mass  flux  of  oxygen 


Fig.  5.  Comparison  of  (a)  diffusive  oxygen  mass  flux  and  (b)  convective  oxygen 
mass  flux  at,  =  0.55  V  (kg  m-2  s_  1 ). 


on  a  Pentium  IV  PC  (2.8  GHz,  2  GB  RAM)  using  Windows  XP 
operating  system. 

3.  Results  and  discussion 

Fig.  3  shows  the  effect  of  overpotential  on  the  gas  velocity 
distributions  in  the  porous  cathode.  Clearly,  the  vectors  direct  the 
gas  flow  from  the  interface  of  the  PEM  and  the  cathode  toward 
to  the  orifice.  In  addition,  a  higher  overpotential  with  a  stronger 
electrochemical  reaction  has  larger  outward  flow  velocities  in  the 
cathode.  The  outward  flow  seems  be  contrary  to  the  common 
sense  that  the  fresh  air  should  flow  into  the  cathode  from  the 
ambient.  Such  transport  phenomena  are  illustrated  as  follows.  In 
the  present  geometric  design,  the  orifice  acts  as  not  only  an  inlet 
of  the  reactants  (O2)  but  also  an  outlet  of  the  products  (H2O).  The 
oxygen  diffuses  into  the  porous  cathode  through  the  orifices  and 
then  is  consumed  continuously  for  keeping  the  oxygen  reduction 
reaction  (Eq.  (1)).  Actually,  the  oxygen  reduction  reaction  will 
increase  the  bulk  mass  of  the  gas  mixture  by  producing  the  water 
vapor.  Thus,  the  gas  mixture  should  flow  outward  to  satisfy  the 
continuity.  It  is  further  seen  from  these  figures  that  significant 


Slice:  Mass  fraction.  w02 


Max:  0.233 

10.232 
0.231 
0.23 
-0.229 
-0.228 
-0.227 
-0.226 
1 0.225 
1 0.224 


Min:  0.223 
.232 

0.22 

0.2 

0.18 

0.16 

0.14 

0.12 

0.1 


Min:  0.0803 

Slice:  Mass  fraction,  w02  Max:  0-232 


Fig.  6.  Effect  of  the  total  overpotential  on  the  mass  fraction  distributions  on 
the  interface  of  electrolyte  membrane  and  catalyst  layer  (z  =  0.000075  mm),  (a) 
rit  =  0.35  V,  (b)  0.55  V  and  (c)  0.70  V. 
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velocity  peaks  occur  around  the  edge  of  the  orifice,  as  denoted  by 
star  symbols  in  each  plot.  This  is  a  contribution  of  the  products 
(water  vapor)  in  the  reactive  layer  beneath  the  current  collector. 

Fig.  4(a  and  b)  show  the  iso-surfaces  of  the  oxygen  mass 
fraction  (coo2)  and  the  water  vapor  mass  fraction  (&>h2o)  in  the 
porous  cathode,  respectively.  The  oxygen  mass  fraction  in  the 
porous  cathode  under  the  orifice  is  almost  constant.  It  decreases 
from  the  orifice  to  the  remote  corner  of  the  domain,  i.e.,  v  =  0, 
y  =  0.  In  contrast,  the  water  vapor  mass  fraction  has  an  inverse 
trend  that  the  remote  corner  has  higher  values  of  cou2o  and  the 
orifice  has  lower  values  of  &>h2o-  In  addition,  both  the  oxygen 
mass  fraction  and  the  water  vapor  mass  fraction  seem  to  be 
unaltered  along  the  z-direction. 

Two  kinds  of  mass  fluxes  in  the  porous  cathode  can  be  decom¬ 
posed  from  the  total  oxygen  mass  fluxes,  i.e.,  diffusive  mass 
flux  and  convective  mass  flux.  As  shown  in  Fig.  5(a),  the  oxy¬ 
gen  diffusive  flux  directs  inward,  which  is  driven  by  the  oxygen 
concentration  gradient.  In  contrast,  the  oxygen  convective  flux 
directs  outward  that  follows  the  gas  flow  direction  that  shows  in 
Fig.  3.  It  is  seen  that  the  magnitude  of  the  diffusive  flux  is  signifi¬ 
cantly  larger  than  that  of  the  convective  flux.  Form  the  integrated 
results  of  the  quarter  of  the  orifice  the  diffusive  mass  flux  for 
the  oxygen  is  about  300  times  of  convective  one.  Therefore, 
the  diffusion  dominates  the  mass  transports  in  the  air-breathing 
cathode. 

Fig.  6(a-c)  show  the  distributions  of  the  oxygen  mass  frac¬ 
tion  (coo2)  at  the  interface  between  the  PEM  and  the  cathode  for 
three  different  total  overpotentials,  i.e.,  rjt  =  0.35, 0.55  and  0.7  V, 
respectively.  It  is  seen  from  these  figures  that  the  projective  area 
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from  the  orifice  has  a  high  oxygen  mass  fraction.  It  decreases 
as  the  oxygen  gradually  diffuses  into  the  cathode  beneath  the 
current  collector.  The  remote  corner  (v  =  0,  y  =  0)  has  a  local 
minimum  mass  fraction  for  oxygen  since  it  has  the  longest  dis¬ 
tance  from  the  orifice  and  thus  cannot  sense  the  fresh  air  easily.  In 
general,  the  oxygen  mass  fraction  decreases  with  increasing  the 
total  overpotential.  At  the  highest  total  overpotential  (Fig.  6(c)), 
the  oxygen  is  depleted  almost. 

Fig.  7(a  and  b)  display  several  slices  cutting  across  the  v 
and  y  planes  showing  the  distributions  of  the  electrolyte-phase 
potential  (<p\)  and  the  catalyst-phase  potential  (0S)  in  the  compu¬ 
tational  module,  respectively.  No  data  of  catalyst  potential  are 
shown  in  the  electrolyte  layer,  since  it  is  free  of  the  catalyst.  It 
is  seen  in  Fig.  7(a),  the  potential  of  the  electrolyte-phase  has 
the  highest  value  on  the  bottom  surface  of  the  module  (z  =  0). 
It  decreases  along  the  z-direction  due  to  the  ohmic  resistance. 
In  addition,  a  significant  drop  of  the  electrolyte-phase  potential 
is  found  near  the  orifice  edge.  It  is  because  a  large  amount  of 
ionic  current  should  traverse  here  to  transfer  the  current  to  the 
catalyst-phase  through  the  electrochemical  reaction.  As  for  the 
potential  catalyst-phase  (Fig.  7(b)),  the  region  under  the  current 
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Fig.  7.  Comparison  of  electrolyte  potential  and  catalyst  potential,  (a)  electrolyte 
potential  (<p\)  and  (b)  catalyst  potential  (</>s)  (unit:  V). 


Fig.  8.  Potential  distributions  along  thickness  direction  (z)  at  x  =  0.75  mm  and 
y  =  0.75  mm:  (a)  electrolyte-phase  and  (b)  catalyst-phase. 
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collector  has  null  potential  almost,  meaning  that  the  porous 
cathode  well  connects  to  the  current  collector  and  the  ohmic 
resistance  is  negligible  almost.  In  addition,  a  longer  distance  of 
through-flow  current  from  the  region  under  the  orifice  to  the 
current  collector  requires  a  higher  electric  potential. 
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Fig.  8  shows  the  potential  distributions  along  the  z-direction 
at  the  module  middle  (x  =  0.750  mm,  v  =  0.75  mm).  Focus  is 
first  placed  on  the  distribution  of  the  electrolyte-phase  poten¬ 
tial  (Fig.  8(a)).  A  linear  drop  of  the  electrolyte-phase  potential 
between  z  =  0  and  0.75  mm  reflects  a  constant  ohmic  resistance 
in  the  PEM.  In  the  porous  cathode,  0.75  mm  x  z  x  1.5  mm,  the 
electrolyte-phase  potential  drops  sharply  first  and  then  declines 
mildly.  The  sharp  drop  of  0 1  is  because  the  through-flow  area 
for  the  ionic  current  is  reduced  when  it  enters  the  porous  cath¬ 
ode  from  the  PEM.  As  the  ionic  current  moves  forward  to  the 
current  collector  (z  =  0. 15  mm),  the  ionic  current  transfers  itself 
gradually  to  the  catalyst  by  the  electrochemical  reaction  and  thus 
eases  off  the  drop  of  the  electrolyte-phase  potential.  In  a  similar 
transport  manner,  the  potential  distribution  of  the  catalyst-phase 
starts  with  a  mild  decline  and  then  drops  significantly  to  zero  at 
the  interface  with  the  current  collector. 

In  the  porous  cathode,  the  difference  between  the  electrolyte- 
phase  potential  and  the  catalyst-phase  potential,  <p\  —  0S,  repre¬ 
sents  the  cathodic  overpotential,  including  the  activation  over¬ 
potential  and  the  concentration  overpotential.  Fig.  9  depicts  the 
local  cathodic  overpotential  distributions  on  three  selected  ele¬ 
vations  of  the  porous  cathode,  i.e.,  z  =  0.075,  0.1  and  0.15  mm, 
respectively.  It  is  seen  that  the  local  cathodic  overpotentials  are 
rather  even  in  throughout  the  cathode.  At  a  fixed  elevation,  the 
cathodic  overpotential  in  the  region  beneath  the  orifice  is  signif¬ 
icant.  In  addition,  the  cathodic  overpotential  is  higher  near  the 
PEM  than  near  the  current  collector.  It  is  because  the  oxygen 
concentration  is  lower  and  electrochemical  reaction  that  requires 
a  larger  overpotential  to  drive  the  transfer  current  between  the 
electrolyte  and  the  catalyst. 
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Fig.  9.  Overpotential  (<j>\  —  <f>s)  distributions  at  three  elevations:  (a) 
z  =  0.075  mm,  (b)  z  =  0. 1  mm  and  (c)  z  =  0.15  mm  (unit:  V). 


Fig.  10.  Ionic  current  density  distributions  and  electronic  current  density  distri¬ 
butions  (unit:  A  cm-2). 
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The  electrochemical  reaction  rate,  represented  by  the  local 
current  density,  is  related  to  both  the  local  overpotential  and  the 
oxygen  concentration.  Fig.  10(a  and  b)  show  the  vector  distri¬ 
butions  of  the  ionic  current  density  and  the  electronic  current 
density  in  the  computational  module,  respectively.  The  ionic 
current  directs  monotonously  from  the  bottom  to  the  top  of  the 
module.  In  contrast,  the  electronic  current  density  distribution 
shows  a  large  variation.  It  directs  toward  the  current  collector 
from  everywhere  of  the  cathode.  Again,  the  current  production 
by  the  electrochemical  reaction  is  higher  beneath  the  orifice  that 
causes  a  high  current  density  in  the  region  around  the  edge  of 
the  orifice. 

Fig.  1 1  shows  the  distributions  of  the  ionic  current  density  and 
the  electronic  current  density  along  the  z-direction  at  the  mid¬ 
dle  of  the  computational  module  (x  =  0.75  mm,  y  =  0.75  mm). 
As  shown  in  Fig.  11(a),  the  ionic  current  density  does  not  alter 
in  the  PEM  but  drops  significantly  in  the  porous  cathode.  The 
extent  of  the  current  density  drop  represents  the  degree  of  the 
electrochemical  reaction.  It  becomes  zero  at  the  interface  with 


Fig.  11.  Current  density  distributions  along  the  z  direction  at  x  =  0,  y  =  0:  (a) 
ionic  current  density  and  (b)  electronic  current  density  (unit:  Am-2). 


Fig.  12.  Polarization  curve  of  the  air-breathing  cathode. 


the  current  collector.  As  for  the  electronic  current  density,  it  tra¬ 
verses  through  the  catalyst  in  the  porous  cathode,  and  increases 
downstream  due  to  the  current  accumulation  that  is  transferred 
from  the  electrolyte. 

Fig.  12  shows  the  polarization  curve  of  the  present  air- 
breathing  cathode.  The  data  shown  in  the  vertical  coordinate 
is  represented  by  the  value  that  subtracts  the  total  overpoten¬ 
tial  from  the  open  circuit  voltage  of  the  fuel  cell.  If  the  anodic 
overpotential  is  negligible,  the  above  data  can  stand  well  for 
the  cell  potential,  Vceii.  It  is  seen  that  cell  potential  drops  sig¬ 
nificantly  at  the  lower  current  density  and  decreases  mildly  as 
the  current  density  increases,  which  is  typical  of  the  PEM  fuel 
cell. 

4.  Concluding  remarks 

A  three-dimensional  model  of  coupled  fluid  flow  field,  mass 
transport  and  electrochemistry  in  an  air-breathing  cathode  of  a 
PEM  fuel  cell  is  presented.  Conservative  equations  describing 
the  co-transports  of  mass,  momentum,  species  and  charge  are 
developed.  A  finite-element  based  CFD  method  is  employed 
to  solve  the  above-coupled  equations  with  proper  account  of 
electrochemical  kinetics.  The  model  has  implemented  a  voltage- 
to-current  algorithm,  coupling  of  the  potential  field  with  the 
species  concentration  field.  A  realistic  spatial  variation  of  elec¬ 
trochemical  kinetics  is  thus  obtained  by  simultaneously  solving 
the  electron  transport  equations  of  the  catalyst  and  electrolyte 
in  the  porous  cathode.  The  effect  of  the  total  overpotential 
on  the  distributions  of  reactant  gas  mass  fraction,  overpoten¬ 
tial  and  current  density  in  the  porous  cathode  are  presented 
and  discussed  in  detail.  They  can  help  in  understanding  of  the 
gas-transport  mechanisms  and  electrochemical  characteristics 
in  an  air-breathing  PEM  fuel  cell,  and  thus  provide  a  solid 
base  for  optimizing  the  geometry  of  the  passive  PEM  fuel  cell 
stack. 
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